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5-H ' We study the effect of impurity scattering in the unitarity limit on the A and B phase superconductivity in 

V^l PrOs4Sbi2. We take the triplet p+h-wave superconducting order parameters and the impurity scattering is 

3 , treated within the standard theory. We find the quasiparticle density of states and thermodynamics are very 

sensitive to the impurity scattering. The impurity scattering dependence of the superconducting transition 
temperature and the superconducting order parameter at T = -R" are very similar to those in d-wave su- 
perconductors. Some of these characteristics will provide a sensitive test of the symmetry of the underlying 
superconducting order parameters. 
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1 Introduction 



Superconductivity in the cubic heavy-fermion (HF) filled skutterudite PrOs4Sbi2 was discovered in 2002 
by Bauer et al|X^, and since that time many experimental and theoretical studies of this compound have 
.^ ' been performed. Vollmer et al O studied the low-temperature specific heat of this compound, finding a 

■r^lj- . double superconducting transition with Tc's of 1.75 and 1.85 K. This material thus appears to have at least 

^D ' two thermodynamic phases, which we term the A and B phases in the following. Vollmer et al also found 

^^ , the value of the specific heat jump AC/C at Tc to be 3, including both transitions. 

Much effort has centered on establishing the phase diagram of PrOs4Sbi2, with several measurements 
S . of the upper critical field Hc2 [3. 4.^- This material is an extreme type 2 superconductor 1 10|, so that in 

the intermediate state between Hd and Hc2 partial flux penetration in the form of vortices appears. All Hc2 
measurements indicate an upper critical field near T=0 of approximately 2 Tesla. Only two measurements 
'^ ' of the A-B phase boundary have been made. Izawa et al 1 6 1 performed a measurement of the magne- 

5 . tothermal conductivity of PrOs4Sbi2 and found that the higher-temperature, higher-field A phase contains 

at least four order parameter nodes along the cubic axes, and that the lower-temperature, lower-field B 
phase contains only two nodes along the [010] axis. Measson et al [7J detected the A-B phase boundary 
through thermodynamic measurements of the specific heat (which undergoes a jump in the transition from 
one phase to another) and initially produced an alternative phase diagram, in which the A-B phase bound- 
Vh ' ary lies much closer to the upper critical field Hc2, and essentially parallels this curve. Recent work by 

Measson, however details several samples which appear to contain only an A phase, so that the phase 
diagram of this material remains highly controversial. Below we present two alternative diagrams based 
upon the work of Measson et al and Izawa et al. 

Much effort has also centered on the determination of the order parameter symmetry of these phases. 
Kotegawa 1 9 1 investigated the nuclear spin-lattice relaxation time via nuclear quadrupolar resonance ex- 
periments, and found evidence for exponentially activated behavior above 0.3Tc, indicating isotropic or 
s-wave pairing. Kotegawa also found a large value of A/Tc of approximately 2.5, significantly exceeding 
the BCS s-wave weak-coupling value of 1.76. Similarly, MacLaughlin et al fTOl studied the pairing state 

pss data will be provided by the publisher 



D. Parker, S. Haas, and K. Maki: Resonance impurity scattering in superconductivity in PrOs4Sbi2 




0.5 



T(K) 



1.5 



.0 0.5 1.0 1.5 
r(K) 



2.0 



Fig. 1 Phase diagrams by Measson et al (left) and Izawa et al. The top hne and data points represent Hc2 of the A 
phase, while the lower one represents the AB phase boundary H*. 



via muon-spin resonance and also found evidence for an isotropic gap. Suderow et al 1111 performed scan- 
ning tunneling microscopy (STM) measurements and asserted that the Fermi surface is gapped over most, 
if not all, quasiparticle momenta. As a directionally locked probe, however, this measurement may have 
difficulty determining certain nodal structures, such as point nodes. 

Other experiments show strong evidence for nodal superconductivity in PrOs4Sbi2. In particular, Chia 
et al fl_2] measured the superfluid density of single crystals and found low-temperature power-law behavior, 
with ps ^ 1 — aT^. This power-law behavior was observed for fields parallel to each of the three crystal 
axes. In addition, the previously mentioned magnetothermal conductivity data of Izawa et al (6| strongly 
suggests the presence of point nodes. 

Several theoretical proposals have been put forth in an attempt to understand the superconductivity in 
PrOs4Sbi2 . Miyake et al 1 1 3 1 proposed triplet p-wave pairing caused by quadrupolar fluctuations, whereas 
Goryo 0| proposed an "anisotropic s-wave" state for the A phase and an "anisotropic sH-id-wave" state 
for the B phase, assuming singlet pairing. Sergienko and Curnoe 1 15 1 proposed states with a mixture of 
line nodes and point nodes. For the purposes of this paper, we will assume the existence of unconventional 
superconductivity in PrOs4Sbi2 with A-phase and B-phase order parameters as given in Ref. |16|, i.e. 
triplet pH-h-wave order parameters. The clean-limit thermodynamics of these order parameters was worked 
out (as in IJSJ) in |,17,|, while the upper critical field Hc2 was worked out in 1 18 1. In both works reasonable 
agreement between theory and experiment was found, suggesting further analysis. 

Some work concerning superconductivity in substituted samples of PrOs4Sbi2 has been performed. 
Frederick et al [201 and Chia et al f2Tl both studied the superconducting order in single crystals of 
Pr(Osi_a;Ru2;)4Sbi2. Frederick found evidence for unconventional superconductivity in pure PrOs4Sbi2, 
but observed exponentially activated behavior in all of the other samples with x > 0.05. Chia found ev- 
idence for a three phase scenario at low dopings (x < 0.2), and confirmed the exponentially activated 
behavior found by Frederick. Rotundu et al 1221 examined the low-temperature susceptibility and specific 
heat of single crystals of Pri_a;Lai_a:Os4Sbi2 and found that the Lanthanum substitution had little effect 
on the superconductivity, with an essentially linear decrease in Tc as x as increased. 

The recent interest in substituted samples suggests that theoretical modeling of the effects of impurity 
scattering is in order We expect the theory outlined below, which assumes strong impurity scattering, 
should be applicable when samples doped towards non-superconducting materials (unlike PrRu4Sbi2 and 
LaOs4Sbi2) are tested. 
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Fig. 2 The quasiparticle density-of-states for the A phase (left) and the B phase are shown, for several impurity 
scattering rates. 



2 Impurity Scattering in PrOs4Sbi2 

In this section we present the results of calculations of the effects of impurity scattering on the supercon- 
ductivity in PrOs4Sbi2. As mentioned previously, we use here the order parameter model described in 
Ref. L16J As we expect the effects of impurity scattering to be quite strong we model the scattering in 
the unitarity limit, in which the phase shift between the incoming and scattered waves is taken to be 7r/2. 
Then following the approach described in Ref. [231 we incorporate the scattering by renormalizing the 
quasiparticle frequency a) as follows: 



zr 



^f^ 



A2/2 



(1) 



Here F is the quasiparticle scattering rate in the normal state and A = A(T). 
1(1 - k^ - ky - kl), and (...) denotes J dQ/An. For the B-phase, / = 1 - 



lo '^^ ■ 



For the A-phase / = 

z"', and (...) denotes 



The quasiparticle density-of-states (DOS) can be obtained once ui is known, from the relation 



N{E) = NaRe < 



V^ 



A2/2 



(2) 



Below is the density-of-states for the A and B phases for a few impurity scattering rates, as indicated. 
As is evident from the plots, even a small amount of impurity scattering creates significant low-energy 
excitations and truncates the peak normally observed at E = A. It is interesting that N(E) is largely 
independent of E for relatively low energy. These DOS should be experimentally accessible via scanning 
tunneling microscopy (STM) measurements. 

Also of interest is the overall dependence of the zero-energy DOS on impurity scattering rates, along 
with the zero-temperature behavior A(F, T = 0) and Tc(F). These quantities are shown below for both 
phases. As stated previously, the zero-energy DOS increases quite rapidly for small impurity concentra- 
tions. Here Tc(r) is found by solving the well-known Abrikosov-Gorkov equation 1241 : 



ln(T,/r,o) 



*( 



,1 



27rr, 



nl) 



(3) 
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Fig. 3 The zero-energy DOS, zero-temperature A(r) and Tc(r) for the A phase (left) and the B phase are shown. 
The zero-energy DOS begins at the lower left corner of each plot. 



Here ^(z) is the digamma function. And Tc, the critical impurity scattering rate where superconductivity 
disappears, is given by 0.8819 Tco for both phases. 

The temperature-dependent order parameter A(r) is determined by solving the BCS weak-coupling 
gap equation: 



Here a;„ and a;„ are related by 



(4) 



Ulri. — t^ri 



v/'^^ + A2/2 



(5) 



Here w„ is the Matsubara frequency ujn — {n + ^)2'kT. In addition, the sum over n must be truncated at 
uj = Eq, the cut-off energy. Plots of the temperature dependent order parameters A a (T) and As (T) are 
presented below. 

Once A(r, T) has been determined, the thermodynamics of the system can be analyzed. We begin with 
the entropy: 



POO 

Ss = -4/ dE N {E){f In, f + {l-f)ln{l-f)) 
Jo 



(6) 



Here / is the Fermi-Dirac distribution (1 + e^^) ^ with /3 ~ l/fc^T and N(E) is the quasiparticle density 
of states. From the entropy the electronic specific heat is determined via 



C -T^ 



(7) 



We show above the specific heat for both phases for all temperatures below Tc for a few impurity scattering 
rates. As expected, the specific heat is T-linear for small T, with the coefficient equal to the zero-energy 
DOS. The specific heat jump AC at Tc decreases rapidly with impurity concentration for both phases. 
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Fig. 4 The temperature dependent order parameters Aa (T) (left) and As (T), for T/Tc = 0, 0.25, 0.5 and 0.75. 



1.5 



H 

t:^ 1 



0.5 



r/r=().25 

r/r=o.5 

r/r=o.75 

Normal state 

r/r=o 



2.5 r 



H 1.5 





— r/r=o.25 




Nonrial state 




____ r/r=o.5 




- - - r/r =0.75 




r/r=o 



o 



0.5 



0.2 0.4 0.6 0.; 
T/T„ 

cO 



"0 0.2 0.4 0.6 0.: 

cO 



Fig. 5 The specific heat for the A phase (left) and B phase is shown, for impurity scattering rates F/Fc = 0, 0.25, 0.5 
and 0.75. 



The superfluid density is calculated as follows; 



,abiT,T) = 2^r5^(-sin2 0- 



21 f\2 



A^l/I 



(c^2+A2|/|2)3/2 



,(r,r) = 27rr^/3cos 



AV\' 



((I;2+A2|/|2)3/2 



(8) 
(9) 



Note that the superfluid density is isotropic for the cubic-symmetry retaining A-phase, but is not for the 
symmetry-breaking B-phase. We show below the superfluid density for both phases for a range of impurity 
concentrations. We observe that while the pure limit superfluid density varies linearly with temperature 
at low temperature for both the A phase and the c-axis B phase results, this dependence changes to a T^ 
dependence when the effect of impurities is considered. 
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Fig. 6 From top left, the superfluid density for the A phase, for the B phase ab plane, and for the B phase c axis. 
Impurity scattering rates are F/Fc = 0, 0.25, 0.5 and 0.75, with lower impurity scattering rates at the top of each plot. 
Note that the top left and bottom plots resemble that of d-wave superconductivity. 



Finally, we consider the nuclear spin lattice relaxation rate T^ ^, which is related to the density-of-states 
N(E) by the following relation: 



/■oo 

{TiT)-^/{TiT)7^T = Tc= / dE N^ {E)sech' 

Jo 



(E/(2T))/2T 



(10) 



Plots of this quantity for both phases, for both the pure and impurity scattering cases are shown. In the 
zero energy limit (TiT)^^ / {TiT)7^^rp becomes A^(0)^. We also note the presence of a small peak in 
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Fig. 7 Plots of the nuclear spin lattice relaxation rate T^ for the A phase (left) and B phase are shown. Impurity 
scattering rates are F/Fc = 0, 0.25, 0.5 and 0.75, with lower impurity scattering rates at the top of each plot. 

both pure case results just below Tc. As will be shown elsewhere [25], this is a general feature of pure- 
case unconventional superconductors obeying BCS weak-coupling theory. It results from the distribution 
of N(E) away from the constant DOS observed in the normal state. Note that even the A-phase, for 
which the density-of-states does not show the usual divergence at E — A, shows a small peak. The B 
phase peak is much larger due to the increased density-of-states near E = A; this divergence, in fact, is 
oc (_E — A)^^/** and is hence a somewhat stronger divergence than the logarithmic divergence found in 
d-wave superconductivity. 
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3 Conclusion 

We have calculated the effect of impurity scattering on the superconductivity in PrOs4Sbi2 within the 
unitary limit. We find that the superconductivity is expected to be quite sensitive to impurity concentration 
within this limit. 
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